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The spin-polarized transport is investigated in a new type of magnetic 

tunnel junction which consists of two ferromagnetic electrodes separated by 
a magnetic barrier and a nonmagnetic metallic spacer. Based on the transfer 
matrix method and the nearly-free-electron-approximation the dependence of 
the tunnel magnetoresistance (TMR) and electron-spin polarization on the 
nonmagnetic layer thickness and the applied bias voltage are studied theo- 
retically. The TMR and spin polarization show an oscillatory behavior as a 
function of the spacer thickness and the bias voltage. The oscillations orig- 
inate from the quantum well states in the spacer, while the existence of the 
magnetic barrier gives rise to a strong spin polarization and high values of 
the TMR. Our results may be useful for the development of spin electronic 
devices based on coherent transport. 
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I. INTRODUCTION 

During recent years, the large magnetoresistance of magnetic tunnel junctions (MTJs) 
has attracted much attention due to the possibility of its application in digital storage 
and magnetic sensor technologies [1,2]. In most of the experiments, the MTJs consist of 
two ferromagnetic metal (FM) electrodes separated by a thin insulator (FM/I/FM). In 
these structures, the electrical resistance is reduced when the magnetization direction of 
the FM layers changes from antiparallel to parallel alignment by an external magnetic field. 
Because the tunnel magnetoresistance (TMR) strongly depends on the spin polarization of 
the magnetic electrodes [3,4], attempts have been made to fabricate junctions with more 
highly spin-polarized ferromagnets [5,6]. 

Recent experiments on the TMR have shown that the insertion of a thin nonmagnetic 
metal (NM) layer between insulating barrier and the FM electrode (FM/I/NM/FM) reduces 
the TMR, when thickness of the NM layer increases [7-10]. This decrease of the TMR can be 
attributed to the decoherence tunneling of electrons from one electrode to the other through 
the NM layer [11]. On the other hand, recent observations by Yuasa et al [12] show clear 
oscillations of the TMR in high-quality NiFe/Al203/Cu/Co junctions in which the Cu/Co 
is a single crystal. This oscillatory behavior has been interpreted in terms of the formation 
of spin-polarized resonant tunneling. 

The effect of nonmagnetic metallic interlayers has been theoretically investigated in 
FM/NM/I/NM/FM structures. Based on the Kubo formalism, Vedyayev et al. [13] studied 
magnetoconductivity as a function of NM layer thickness at low bias voltages. By using the 
transfer matrix approach, Wilczynski and Barnas [14] calculated the thickness and voltage 
dependence of TMR. These calculations showed that the presence of thin NM spacers can 
lead to the formation of quantum well states that lead to oscillations of magnetoconductivity 
and TMR. When the thickness of only one of the NM layer varies, and also in the case of 
an asymmetric structure, such as FM/NM/I/FM [15], the sign of magnetoconductivity and 
TMR is predicted to oscillate as a function of the NM layer thickness. 



TMR has also been investigated both experimentally [16,17] and theoretically [18-21] in 
spin filter tunnel junctions in which a ferromagnetic semiconductor (FMS) is applied as a 
magnetic barrier. The large spin polarization achievable using magnetic barriers [22] makes 
spin filtering a nearly ideal method for spin injection into semiconductors, enabling novel 
spintronic devices [23]. When an FMS layer which acts as a spin filter is used in tunnel 
structures, due to the spin splitting of its conduction band, tunneling electrons see a spin- 
dependent barrier height. Thus one spin channel will have a larger tunneling probability 
than the other, and a highly spin-polarized current may result. 

The purpose of the present work is to investigate theoretically the effect of a magnetic 
barrier on spin transport in the presence of a nonmagnetic spacer layer. Using the transfer 
matrix method and the nearly-free-electron approximation we will study the variation of 
TMR and spin polarization of tunneling electrons in terms of the spacer thickness and bias 
voltage at T=0 K. We assume that the electron wave vector parallel to the interfaces and 
the electron spin are conserved in the tunneling process through the whole system. 

The paper is organized as follows. In section 2, the model and the analytic formulae of 
the spin polarization and TMR for FM/FMS/NM/FM structure are presented. In section 
3, numerical results obtained for a typical tunnel junction are discussed. The results of this 
work are summarized in section 4. 

II. DESCRIPTION OF MODEL AND FORMALISM 

We consider two semi-infinite FM electrodes separated by a ferromagnetic barrier which 
acts as a spin filter, and a nonmagnetic metallic layer as a quantum well, in the presence 
of DC applied bias Va as shown in Fig. 1. For simplicity, we assume that the two FM 
electrodes are made of the same metal and all of the interfaces are flat. In a nearly-free- 
electron approximation of spin-polarized conduction electrons, the longitudinal part of the 
effective one-electron Hamiltonian can be written as 



H^ 



2m* dx'^ 



+ U,j{x)-\iya + V^ 



(1) 



where m* (j=l-4) is the electron effective mass in the jth layer, and 



U,{x) 



(2) 



0, X < , 

^FL + - eKa;/tbar, < X < tbar , 

— eVa, tbar < X < tbar + tNM , 

— eVa, X > tbar + tNM , 

where E-ph is the Fermi energy in the left electrode. — hj ■ a is the internal exchange energy 
where \\j is the molecular field in the jth FM electrode and a is the conventional Pauli spin 
operator. The last term in Eq. (1) is a spin-dependent potential and denotes the s — f 
exchange coupling between the spin of tunneling electrons and the localized / spins in the 
magnetic barrier. This term, within the mean-field approximation, is proportional to the 
thermal average of the / spins, (5*2), and can be written as V^ = —aI{Sz). Here, a = ±1 
which corresponds to a =t, i, and / is the s — f exchange constant in the magnetic barrier. 
The Schrodinger equation for a biased barrier layer can be simplified by a coordinate 
transformation whose solution is the linear combination of the Airy function Ai[p(x)] and its 
complement Bi[p(x)] [24]. Considering all four regions of the FM/FMS/NM/FM junction 
shown in Fig. 1, the eigenf unctions of the Hamiltonian (1) with eigenvalue E^ have the 
following forms: 

A2aP^l[pa{x)] + fi2aBi[p^(x)], < X < tbar , 

Ag^e^'^^^ + BsaC-''''''' , tbar <X < tbar + tNM , 

V 

where. 
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(5) 



h^ = J2ml{E^ + eVa + rihQa)/Ti , 



(6) 



are the electron wave vectors along the x axis. Here, ho = \hj\ and r] = +1(— 1) for parallel 
(antiparallel) alignment of the magnetizations. The coefficients Aj„ and Bj„ are constants 
to be determined from the boundary conditions, while 
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Pa{x) = J + f^'^ ^ 
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Although the transverse momentum k|| is omitted from the above notations, the summation 
over k|| is carried out in our calculations. 

Upon applying the boundary conditions such that the wave functions and their first 
derivatives are matched at each interface point Xj , i.e. ipj^rrixj) = ipj+i^a{xj) and 
{m*)^^[diljj^cr{xj)/dx] = {m*^i)^^[diljj+i^cr{xj) /dx], we obtain a matrix formula that connects 
the coefficients Ai^ and Bi^ with the coefficients A^^ and B^^ as follows: 



where 



4cr 



Mtotal — -, 

ha 



Axa 






M 



total 



Aa, 

B4a 



(10) 






A ml 



Ai[p^(a; = 0)] Bi[p^(a; = 
Ai'[p^(a; = 0)] Bi'[p<^(a; = 

Ai[p^(x = tbar)] Bi[p^(x = tbar)] 

iiAi'[/'<^(^ = ^bar)] ^;^Bi'[p^(a; = tbar)] 



-1 



X 



X 



cos(A;3tNM) -^sin(A;3tNM) 




iki^^ m\ 


^ sin(A;3tNM) cos(A;3tNM) 




ik^cF —ml 


g — ifc4CT(tbar+tNM) Q 


-1 


g g«fc4o-{tbar+tNM'\ 





-1 



(11) 



Since there is no reflection in region 4, the coefficient -840- in Eq. (3) is zero and the 
transmission coefficient (TC) of the spin a electron, which is defined as the ratio of the 
transmitted flux to the incident flux can be written as 

1 



T.(E.,K)-^''^'^^ 



^total 



(12) 



ki„ml 

where M}^^^i is the left-upper element of the matrix Mtotai defined in Eq. (11). At T=0 K, 
the spin-dependent current density for the magnetic tunnel junction at a given applied bias 
Va within the nearly- free-electron model is given by the formula [25]: 
em* 



JaiVa) 



A-K^h^ 



eVa / T^{E,, Va)dE, + / (Efl - E^)T„{E,, Va)dE, 

JE^ JEwT-eVr, 



, (13) 



where E^ is the lowest possible energy that will allow transmission and is given by Eq = 
max{— /iq, —eVa — i]hQ} for spin-up and Eq = Hq for spin-down electrons. It is clear that the 
tunnel current is modulated by the magnetic configurations of the both FM electrodes. 
The degree of spin polarization for the tunnel current is defined by 



J^-Ji 



(14) 



J^ + h 
where J| (J|) is the spin- up (spin-down) current density. For the present structure, we 

obtain this quantity when the magnetizations of two FM electrodes and the FMS layer 

(magnetic barrier) are in parallel alignment. For calculating the TMR, both spin currents 

in the parallel and antiparallel alignments are calculated. In this case, the TMR can be 

defined as 



( jf + jf) - {jf + jr 



TMR= '-'--'(p ;'ap ' ' ■ (15) 



where J^ ^ corresponds to the current density in the parallel (antiparallel) alignments of 
the magnetizations in the FM electrodes, for a spin- up (t) or spin-down (|) electron. 

III. NUMERICAL RESULTS AND DISCUSSION 

Numerical calculations have been carried out to investigate the effect of a magnetic 
barrier and an NM metallic layer on the spin transport in the Fe/EuS/Au/Fe structure as 
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a typical MTJ. We have chosen Fe and EuS because they have cubic structures and the 
lattice mismatch is very small [26]. The parameters E-p and Hq for Fe electrodes are taken 
corresponding to A;f|=1.09 A^^ and kpi=0A2 A~^ (for itinerant d electrons) [27]. For the 
EuS layer we use S=7/2, 1=0.1 eV [28], 0=1.94 eV [21] as a symmetric barrier height, and 
tbar=l-3788 nm which corresponds to four monolayers of EuS(lll) [29]. The Fermi energy 
in the Au layer is -Ef.nm = 5.51 eV [30]. In practice, the effective masses of electrons may 
differ from that of free electron, but here for simplicity, we assume all electrons have the 
same mass, rrie as free electrons. 

In our considered system, the magnetization direction of the left Fe electrode and the EuS 
magnetic barrier stays fixed, because they are coupled through the interface. Inserting the 
NM spacer layer between the FMS layer and the right FM electrode decreases the exchange 
coupling between the two FM electrodes. Therefore the right Fe electrode is free and may 
be switched back and forth by an external magnetic field (see Fig. 1). In the following, we 
show the numerical results at T=0 K, so in the magnetic barrier we have (Sz) = S. 

The TC for tunneling electrons at the Fermi energy, and the TMR in the tunnel junction 
as a function of the NM layer thickness are shown in Figs. 2(a) and 2, respectively. It can 
be seen that both the TC and the TMR oscillate with the NM layer thickness due to the 
dependence of these quantities on the incident electron's energy and the NM layer thickness. 
The physics behind this behavior is easy to understand. When the energy of the incident 
electrons is approximately equal to the energy of a quasibound state in the quantum well, a 
resonance condition is fulfilled and the TC of the electrons through the structure strongly 
increases. On the other hand, the position of the quantum well levels, formed in the NM 
layer, strongly depends on the well thickness. Therefore, with continuous variation of the 
tNM, the position of the resonant states varies and this leads to the oscillations of the TC 
and the TMR. From Fig. 2, it is clear that where the transmission is high, the TMR is low. 
From the continuous variation of tNM, one can see that both the TC and the TMR oscillate 
with a period of 0.26 nm which is in excellent agreement with Xpf^, where Xp is the Fermi 
wavelength in the Au layer and is equal to 0.522 nm. The results show that, for the spin-up 
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electrons, the TC is very higher than the spin-down ones both in the parallel and antiparallel 
alignment of magnetizations. The reason can be explained by the spin filtering effect of the 
FMS barrier. As shown in Fig. 1, due to the spin splitting of the conduction band in the 
FMS layer, the barrier height for spin-up electrons is lowered and the tunneling probability 
greatly increases, while for the spin-down electrons the barrier height is raised and hence, 
the tunneling probability for these electrons reduces. Thus, a large spin-polarized current is 
produced. 

In Fig. 2(b) we have also shown the TMR in terms of thickness, when tNM changes 
in monolayer steps. At the monolayer thickness tML,NM=0.2355 nm, which corresponds to 
interlayer distance of Au(lll), the TMR first increases in the first three monolayer and after 
that the TMR becomes less and starts to oscillate with a period of approximately 2.5 nm. 
In contrast, for tML,NM=0.27 nm, the TMR initially decreases and the period of oscillations 
is approximately 7 nm. 

The voltage dependence of the TMR for several thicknesses is shown in Fig. 3, when 
^ML,NM=0.2355 nm. We can see that, due to the FMS layer, the TMR is asymmetric, as 
a tunnel junction without a nonmagnetic layer is used. For each fixed voltage, the TMR 
oscillates, as is shown in Fig. 2(b) for \4=50 mV. The results show that more than 250% 
TMR can be obtained in reverse bias. By increasing the number of nonmagnetic monolayers, 
the TMR shows an oscillatory behavior in the positive voltages. Under a forward bias, 
electrons at the Fermi energy of the left electrode, tunnel into the empty states of the 
quantum well, which have been created due to the applied bias. Therefore the applied bias 
has a strong influence on the TMR and the oscillatory behavior will appear. On the other 
hand, under a reverse bias, electrons at the Fermi energy of the quantum well tunnel into 
the empty states of the left electrode. In this case, variations of the applied bias cannot 
strongly affect the TMR. 

In order to further see the effect of the FMS layer on the tunnel currents, we have shown 
in Fig. 4 the spin polarization of tunneling electrons as a function of the bias voltage, for 
^ML,NM=0.2355 nm. It is obvious that the spin polarization is symmetric for tML,NM=0, and 
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becomes strongly asymmetric when the number of monolayers increases. The results show 
that the spin polarization can reach 99%, which is evidence of spin filtering effect in the 
FMS layer. 

The experimental measurements show that, if one uses a nonmagnetic insulator instead of 
a magnetic tunnel barrier, very low values for electron-spin polarization and TMR will result 
[7,12]. The oscillatory behavior of TMR in terms of the NM layer thickness, is qualitatively 
in agreement with the experiments of Yuasa et al [12], and shows that the coherent tunneling 
of electrons through the quantum well states is the main reason of this behavior. 

The present results can be qualitatively compared with those already known in the 
relevant literature, when a nonmagnetic barrier is used. In [13,14], the authors showed 
that in the FM/NM/I/NM/FM systems, the TMR in two different cases has an oscillatory 
behavior with NM layers: the case where the thickness of both NM layers are equally varied 
(symmetric case), and the case where only one of the NM layers is increased (asymmetric 
case). In the symmetric case, the TMR is always positive, but in the asymmetric case the 
TMR oscillates from positive to negative as a function of the NM layer thickness. The 
situation of the asymmetric case is similar to our structure in which only one NM layer has 
been used [15]. Thus, as it is clear from Figs. 2(b) and 3, our results are in qualitative 
agreement with the asymmetric case in FM/NM/I/NM/FM structures. 

Finally, we discuss the effect of interface scattering on spin-injection from the FM into 
NM spacer [31,32]. For this purpose, we treated the magnetic barrier as a (5-function U(j5{x), 
with U(j describing the barrier strength parameter for spin a electron. In this case, the pa- 
rameter U„ is proportional to both height and width of the barrier. The numerical results 
showed that with increasing the f/o-, the TC reduces, and hence the spin-dependent cur- 
rents decay exponentially, which is due to the interface resistance. On the other hand, the 
spin polarization and TMR increase with f/o-, and trend gradually towards constant values, 
because in these values of the barrier parameter, only the FM electrodes have a dominant 
influence on the spin polarization and TMR. 



IV. CONCLUDING REMARKS 

The magnetotransport behavior in an FM/FMS/NM/FM tunnel junctions was investi- 
gated theoretically, using the transfer matrix method and the nearly-free-electron approxi- 
mation. The numerical results indicate that in Fe/EuS/Au/Fe structure, the TMR has an 
oscillatory behavior as a function of the NM layer thickness. By increasing the thickness 
of the NM layer, both the TMR and electron-spin polarization, oscillate under the forward 
biases. Due to the presence of a magnetic barrier, there exist more than 250% TMR and 
98% spin polarization in the tunnel currents; this can be obtained by adjusting the thickness 
and the applied bias. 

Our study restricted to the zero temperature limit, because the ferromagnetic transition 
temperature in most of the FMS layers is much lower than room temperature [29]. Low 
temperature spin transport is, however, important as a testing ground for novel ideas and 
concepts in future spin electronic devices based on coherent transport [12,33]. 
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FIG. 1. Spin-dependent potential profile for an FM/FMS/NM/FM magnetic tunnel junction 
under forward bias Va- In the FMS layer, the dashed line represents the bottom of the conduction 
band at T > Tc and the thin arrows indicate the bottom of the conduction band for spin-up and 
spin-down electrons at T < Tc- The zero of energy is taken at the middle of the bottoms for the 
majority-spin band and minority-spin one in the left FM electrode. 
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FIG. 2. (a) Spin-dependent TC at Fermi energy in the parallel (P) and antiparallel (AP) con- 
figurations, and (b) dependence of the TMR, as a function of Au layer thickness. 
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FIG. 3. Dependence of the TMR on the apphed voltage for different NM multilayers. 
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FIG. 4. Dependence of the spin polarization on the applied voltage for different NM multilayers. 
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